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ABSTRACT 

It has been demonstrated that the inclusion of baryonic physics can alter the dark matter densities 
in the centers of low-mass galaxies, making the central dark matter slope more shallow than predicted 
in pure cold dark matter simulations. This flattening of th e dark matter profil e can occur in the 
most luminous subhalos around Milky Way-mass galaxies. I Zolotov et al.l (|2012t ) have suggested a 
correction to be applied to the central masses of dark matter-only satellites in order to mimic the 
affect of (1) the flattening of the dark matter cusp due to supernova feedback in luminous satellites, 
and (2) enhanced tidal stripping due to the presence of a baryonic disk. In this paper, we apply this 
correction to the z =0 subhalo masses from the high resolution, dark matter-only Via Lactea II (VL2) 
simulation, and find that the number of massive subhalos is dramatically reduced. After adopting a 
stellar mass to halo mass relationship for the VL2 halos, and identifying subhalos that are (1) likely 
to be destroyed by stripping and (2) likely to have star formation suppressed by photo-heating, we 
find that the number of massive, luminous satellites around a Milky Way-mass galaxy is in agreement 
with the number of observed satellites around the Milky Way or M31. We conclude that baryonic 
processes have the potential to solve the missing satellites problem. 



1. INTRODUCTION 

The cosmological paradigm based on cold dark mat- 
ter (CDM) and dark energy (A) has been wildly suc- 
cessful in describing the observed evolution and large 
scale structure of our Universe. At small scales, how- 
ever, a number of observations seem to be at odds with 
the predictions of ACDM cosmolog y. In particular, over 
a de cade ago it was poi nted out by iMoore et al.1 ()1999l ) 
and iKlvpin et al.l (|1999f ) that the number of high mass 
subhalos predicted by high resolution CDM simulations 
exceeds the observed number of luminous satellites of 
the Milky Way (MW) by at least an order of magnitude. 
This has become known as the "missing satellites prob- 
lem" (MSP). Although this problem has been mitigated 
to a degree by the dis covery of a number of additional 
faint satellite galaxies dWillman et al.l 120051 Ihwin et al 
Liu et aij|2q08t iMartin et al .1120081: iSimon fc Geha 
GrebefeOOOHvan den Berghll2000t (Belokurov et al 
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Photoionization is expected to suppress star form ation 
in halos below -TO 9 M Q (jOkamoto et all 120081) . In 
more massive halos where gas is retained and star for- 
mation can begin, fur ther suppress i on is expected from 
supernova feedback (IDekel fc Silkl [l98l IBenson et al.l 
[2001 iDekel fc Wool [20031: IGovernato et al l I2007D . It is 
also important to consider the number of faint satel- 
lites that re main unobserved due to observational in- 
completeness |Willman et a l. 2004; S imon fc Gehal[2007l: 
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Tollerud et alJ 120081: iWalsh et al.l l2009t iKoposov et al.1 



2009t IRashkov et all 120121) . Taken together, it is pos 



remains a considerable discrepancy between the number 
of observed MW satellites and the number predicted in 
CDM simulations. 

Efforts to resolve this issue have fallen into two broad 
categories. First, there are proposals in which the 
star formation rate in satellite galaxies is suppressed, 
leading to large numbers of low mass subhalos which 
are simply unobservable. Possible means for such sup- 
pression include photoevaporation r esulting from ioniz- 
ing radiation, e.g., at reionization ([Quinn "eFHI [T996I: 
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sible that observational incompleteness combined with 
suppression of star formation can bring the number of 
luminous satellites in line with the number of predicted 
satellites around a Milky Way-mass galaxy. 

However, while this combination of effects may 
explain the faint or low-m ass regime of the MSP, 
iBoylan-Kolchin et al.l (|2011l ) recently pointed out that 
it still fails at the massive end. Dubbed the "too big 
to fail" (TBTF) problem, this aspect of the MSP arises 
because the most massive subhalos in ultra-high resolu- 
tion djixl£j2ji : ttaM3n galax- 
ies (jDiemand et al.l 120081 iSpringel et all 120081 ) are too 
dense to host the obs e rved satellites of the MW (see also 
IWolf fc Bullocg[20Tl iHavashi fc Chiball2012l ). The sim- 
ulations always contain a population of subhalos (6-22, 
varying within the errors of the MWs measured mass) 
that are more massive than any of the dwarf sphero idals 
observed in the MW ijBovlan-Kolchin et al.ll2012[) . In 
other words, while the abundance of lower luminosity 
satellites may be made consistent with CDM predictions, 
the MSP remains a puzzle because simulations predict 
too many massive satellites. 
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A second class of possible solutions to the MSP con- 
siders departures from the standard assumptions of cold 
and collisionlcss dark matter. Whereas typical particle 
dark matter models predict that WIMPs will form ha- 
los with masses as small as 10 -6 Mq or so (depending 
on the temperature at which the WIMPs undergo ki- 
netic decoupling with the cosmic neutrino background), 
the formation of small scale structure can be strongly 
suppressed if the dark matter particles are not en- 
tirely cold rtCohn et all 120001: iMaccio fc Fontanotl l2010t 
lLovell et al.l 1201 2f T! Dark matter in the form of sterile 
neutrinos with masses of ~ 1-10 keV have received at - 
tention within this context (|Dodelson fc Widrowl 119941 ) . 
although many other warm dark matter (WDM) scenar- 
ios could potential ly accomodate a similar suppression of 
small scale power (jPagels & Primack 19821: iHooper et all 
[2001 . WDM primarily addresses the low-mass end of the 
MSP by suppressing the formation of small scale struc- 
ture. It can also help at the high mass end, however, 
because the delayed halo collapse times in a WDM cos- 
mology result in l ower concentration s and hence reduced 
central densities (jLovell et al.H2012T) . 

Self-interacti ng da rk matter (SIDM) models 
(ICarlson et all fl99l ISpergel fc Steinhardti 1200(1 
lLoeb fc Weinerl 120111) provide another mechanism to 
achieve the same effect. In this case, the interactions 
prevent the formation of the steep central density cusps 
that are the hallmark of CDM halos. Density profiles in 
SIDM instead exhibit a c entral core (jVogelsberger et all 
120121 IRocha et all I2012D . which helps to address the 
MSP in two ways: cored halos are more susceptible 
to tidal disruption, potentially removing many of the 
excess low mass halos; secondly, in the surviving halos, 
the core reduce s the central densities, alleviating the 
TBTF problem dVogelsberg cr et al.H2012i ). Furthermore, 
there is considerable observational evidence for the 
existence of dark matter cores in the centers of low 
surface brightness galaxies (iKuzio de Narav et al.l 120081 : 
Ide Blokll2010T: lOh et all [2 01 lft and in at least two MW 
dwarf spheroidal satellites (j Walker fc Peharrubial 120111 ). 

In this work we consider a new type of solution to the 
MSP in which the shape of satellite dark matter profiles 
are altered, but as a result of baryonic physics, rather 
than through modifications in the particle physics sector. 
This model incorporates suppression of star formation 
from photoionization and supernova feedback, but addi- 
tionally considers the tidal effects due to the presence of 
a baryonic disk, which is not found in dark matter-only 
simulations. The loss of gas in satellite halos through 
tidal stripp ing has been propose d as a means to limit star 
formation dStrigari et al.l [20071 ). However, tidal strip- 
ping in a cuspy halo by itself does not reduce the central 
densities of the most massive subhalos enough to bring 
them into agree ment with the observed kinematics of the 
MW satellites ([Read et alj|2006at iBovlan-Kolchin et all 
|2012[). An additional modification to the central densi- 
ties of the most massive satellites appears to be required. 

It has become broadly accepted in recent years 
that baryonic processes can alter the distributions 
of d ark matter withi n halo s, potentially s t eepen- 
ing dBlumenthal et all ITiM iGnedin fc Zhaol [20021: 
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depending on the mass of the halo in question, and on 
the strength of various feedback mechanisms. Recently, 
the effects of baryons on the dark matter halo profiles of 
lumin ous satellite galaxies was studied by iZolotov et all 
(|2012ft . These authors found that satellites with 
Mstar <: 10 7 -M© tend to develop "cored" density distri- 
butions through supernova feedback prior to infall (here, 
cored refers to any inner density slope, 7, shallower than 
-1, becau se 7 becomes flatter w ith increasing stellar 
mass, see IGovernato et al.l I2012T ). The flatter density 
profiles make them more vuln erable to tidal effects 
of the baryonic d i sk (see also iTavlor fc Babul! 120011: 
Read et all l2006al: IChoi et al.l [20091: iWetzel fc White! 
20101: iRomano-Di'az et alll2010al iD'Onghia et alll2010a 
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Penarrubia et all 120101 ) . IBrooks fe Zolotovl Jm% 
showed that the presence of the disk increases tidal 
stripping for galaxies of all masses compared to the dark 
matter-only case, but most strongly reduces the central 
densities of cored satellites (see also iStoehr et all 120021 : 
lHavashi et al.l 120031: iPeharrubia et al.l I2010T ) . which are 
the most luminous. In other words, baryonic physics 
significantly reduc es the observable mass of satellites. 
IBrooks fc Zolotovl (120121 ) demonstrated that these 
combined baryonic effects produce a z ~ satellite dis- 
tribution with luminosities and kinematics comparable 
to the MW and M31. 

In t his work, we apply the results of IZolotov et all 
(|2012L hereafter Z12) to the subhalo populations found 
in the DM-only Via Lactea II simulation. We find that 
after taking into account baryonic effects, the number of 
surviving massive subhalos is strongly reduced. Adopt- 
ing a stellar mass to halo mass relation, we show that 
the number of luminous satellites is also compatible with 
observations. This suggests that a proper accounting of 
baryonic effects reduces the predicted number of massive, 
luminous subhalos that should survive in CDM, bringing 
the observations and theory into agreement without the 
need of invoking warm dark matter or dark matter with 
exotic properties or interactions. 

2. SATELLITE MASSES 

We make use of the publicly available subhalo cata- 
log^ from the Via La ctea II (VL2) DM-only simulation 
(jDiemand et al. 2008}), supplemented by additional infor- 
mation extracted from the simulation outputs, consisting 
of (1) mass enclosed within 1 kpc at z = 0, (2) infall times 
and the maximum value in their rotation c urves, ti„| M , at 
infall , and (3) their full orbital information (jKuhlen et all 
12012ft . including number and distance of apo- and peri- 
center passages. T he VL2 halo was run using a WMAP 
year 3 cosmology (jSpergel et al.l [20071) . The halo was 
selected to have no low z major mergers, thought to be 
similar to the merging history of the MW. The z — halo 
has roughly 450 million particles with masses of 4100 M Q 
within its viriaQ radius of 402 kpc at z = 0, yielding a 
halo mass of 1.9 xlO 12 M . 

Starting from all subhalos within i2 v ; r at z = 0, we 
restrict our study to only those subhalos with i> pea k 

5 http:/ /www. ucolick.org/~diemand/vl/ 

6 Here we adopt a virial overdensity of p/po = 200. 



No Massive Missing Satellites 



3 




3 10 30 100 



v [km s 1 ] 

Fig. 1. — Cumulative velocity functions for VL2 halos that sur- 
vive to z = 0. The ^peak-function is plotted with a solid black line 

for ripeak > 20 km s . The corresponding evolved u max -functions 
are also shown at infall (red, short dashed line) and at z = (dot- 
ted line). For comparison, we also plot i>ik pc at 2: = (blue, long 
dashed line). The ^ pea k an d u max -functions found in the full VL2 
catalog (no cut on f pea k) are shown as continuations in light grey. 

(the largest value of v max over the subhalo's entire his- 
tory) > 20 km s" 1 . It is assumed that all halos with 

Wpoak < 20 km s _1 are st rongly affected by UV heating, 
unbinding their gas (e.g.. IQkamoto et"aT]|2008l) . This ei- 
ther leaves them entirely devoid of stars ( "dark" ) , or very 
faint if they were able to form a few stars prior to reion- 
ization (e.g.. iBovill fc Ricott i 2009; Salvadori fc Ferraral 
[20091 iLi et all 120101: IWolf et al.l I2010f ) . and prone to be 
missed by current surveys. Either way, they are well be- 
low the mass and luminosity of the satellites contributing 
to the massive end of the MSP, and excluding them has 
no impact on our conclusions. In essence, our cut limits 
us to subhalos with masses greater than a few x 10 s M© 
at accretion. This cut yields a sample of 125 subhalos. 

The v max -function of our final subhalo sample (i.e., 
those halos that survive to z = and had i> pea k > 
20 km s _1 ) is shown in Figure 1, for various times. The 
solid line shows the cumulative u pea k-function for the ha- 
los in our sample, while the short-dashed curve shows 
Wmax at infall. We note that w poa k and i> max at infall, 
Winfaii, can be slightly different for these subhalos. The 
majority of the VL2 subhalos in our sample reach u pca k 
at high z and then grow very little, with an overall slight 
decrease (< 10%) in v max . A few subhalos undergo a 
substantial reduction in w max between high z and infall, 
presumably due to encounters with other halos that strip 
their mass. The resulting velocity function at z = is 
shown in Figure 1 as the dotted line. The circular veloc- 
ity at 1 kpc, uik P c, is shown for comparison (long-dashed 
line). In the following, we adopt t>ik pc to examine the ef- 
fect of baryons on the central mass distribution of satel- 
lite galaxies. 

Z12 proposed a correction to the «ik P c values of z = 
DM-only subhalos to account for the missing baryonic 
physics that lowers the central masses of luminous sub- 



halos, 

A(wik P c) = 0.2« infall - 0.26 km s _1 . (1) 

This correction is valid for DM-only halos with 
20 km s™ 1 < Ujnfall < 50 km s _1 . Extrapolating the 
correction to subhalos with i>i n fail < 20 km s - will 
yield unphysical (negative) results for those halos. At 
Vinf a ii > 50 km s _1 , subhalos are Magellanic-like and gas- 
rich at accretion, possibly including an additional ef- 
fect of adiabatic contraction that is not accounted for in 
the correction. There are 5 massive VL2 subhalos with 
Winf a ii > 50 km s" 1 , for which the Z12 correction is not 
applied. 

For subhalos with 20 km s _1 < Vmfaii < 30 km s^ 1 , the 
Z12 correction is designed to account for a reduction in 
central mass due to (1) loss of gas, either due to heating 
before infall or stripping after infall, and (2) the tidal 
effect of the baryonic disk (which does not exist in DM- 
only runs) . These processes should act on a subhalo even 
if it is too low mass to retain gas at reionization, or mas- 
sive enough to be luminous. These processes should also 
occur irrespective of whether the subhalo has a cuspy or 
a cored DM density profile. Yet neither are typically ac- 
counted for in DM-only simulations. In subhalos with 
^infall > 30 km s _1 , Z12 found that enough star forma- 
tion takes place to significantly flatten the DM density 
profiles prior to infall. Hence, at f infall > 30 km s _1 , the 
Z12 correction accounts for (1) and (2) as in the lower 
mass case, but also for an additional reduction in the 
central masses of the satellites and the enhanced tidal 
stripping that occurs as the central density profiles be- 
come more shallow. 

Most of the MW and M31 dwarf spheroidal galaxies 
(dSphs) have central v c values < 20 km s -1 (measured at 
the h alf light radii, which are all < 1 kpc, IMcConnachiel 
120121 ). At z = 0, the VL2 host halo contains 28 satel- 
lites with vik P c > 20 km s 1 , grossly inconsistent with 
the observational results. After applying the correction 
of Z12, there are only 5 satellites with i>ik pc > 20 km s _1 . 
This is a first indication that baryonic physics, at least 
as implemented in the simulations of Z12, appears to be 
quite capable of reducing the mismatch between the cen- 
tral densities of the most massive subhalos in DM-only 
simulations and the observed kinematics in the classical 
MW dwarf satellite galaxies. 

3. LUMINOUS SATELLITES AT REDSHIFT 

It is clear that the baryonic correction proposed by 
Z12 dramatically reduces the number of massive halos 
expected in a DM-only MW-mass run (e.g., from 28 to 
5 for satellites with vik pc > 20km s _1 in VL2). How- 
ever, it is less clear if it reduces the number of luminous 
subhalos. That is, does the correction take 20 Fornax- 
like satellites that initially have wik pc > 20 km s" 1 at 
z = and simply shift them to «ik pc < 20 km s _1 ? If so, 
then there may be fewer massive halos, but the overall 
number of luminous satellites would still be much larger 
than observed in the MW or M31. In this section we 
explore whether both the mass and luminosity of Milky 
Way-mass galaxies can be reproduced by accounting for 
baryonic effects. 
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3.1. Destruction Rates 

We must first consider whether the VL2 halos should 
have all survived to z = 0, or whether it is likely that 
baryonic physics would have led to their destruction. 
iPefiarrubia et all (I2010T ) examined the effect of a bary- 
onic disk on both cored and cuspy subhalos. They 
showed that tides from the disk are dominant in the 
center of cored subhalos. leading to much more mass 
loss than DM-only runs predict. Even cuspy subhalos 
will undergo more mass loss when a baryonic disk is 
present if they are on highly elliptical orbits. Impor- 
tantly, cuspy halos tend to survive even after substantial 
(99.99%) mass loss, while cored halos can be completely 
disrupted. This implies that there could be a large num- 
ber of satellites that survive in the DM-only VL2 run 
that would not survive in a baryonic run. 

In what follows, wc assume that a VL2 subhalo that 
has lost a certain fraction of its mass since infall should 
have been fully disrupted if baryons had been included, 
and should no longer appear as a bound luminous satel- 
lite at z = 0. To set the l i miting fraction, we refer to 
figure 2 in IPefiarrubia et all (|2010l ) . where tidal effects of 
a disk on cored dwarf galaxies are compared to those 
on cuspy halos with no disk present. Adopting their 
"mixed" model, in which a cored (7 = 0) subhalo evolves 
within a parent halo with a cuspy (7 = —1) density pro- 
file, we find that cored subhalos that have pericenters 
< 20 kpc undergo 99.9% mass loss, which we consider 
to be comple tely disrupted. Without the disk presence, 
the results in IPefiarrubia et all (|2010l) show that a cuspy 
subhalo with a pericenter of 20 kpc would experience 
90% mass loss. Hence, we assume that any VL2 subha- 
los with pericenters < 20 kpc should experience strong 
disk tides that need to be accounted for. We conclude 
that if a VL2 DM-only satellite has lost more than 90% 
of its mass after infall, and has had pericentric passages 
that take it within 20 kpc of the parent halo's center, the 
same satellite in a baryonic run is extremely likely to be 
fully disrupted by tides. 

To be conservative, we apply this limit only to subhalos 
that h ave Ui n f a ii > 30 km s" 1 . Z12 and iGovernato et al.l 
(|2012f ) demonstrated that cored DM density profiles ex- 
ist in satellites with a stellar mass of more than 10 7 M , 
corresponding to halos with Uinfaii > 30 km s _1 . While it 
is possible that halos at lower masses have cores (but at 
smaller radii that are below the resolution limit of current 
cosmological simulations), wc assume they may remain 
cuspy, and 90% mass loss may then not be enough to 
fully destroy them. Instead, for the remainder of the 
h alos with i>j n f a n < 30 k m s _1 , we adopt the results 
of I Wetzel fc White! pOlOh and assume that a halo must 
have lost 97% of its mass since infall to be fully stripped, 
and no longer appear as a bound luminous satellite at 
2 = 0. We verified that the majority of the halos that 
have lost 97% of their mass have tidal radii less than 1 
kpc, suggesting that their inner luminous regions should 
indeed be stri pped. Adopting 97% i s conservative, as 
both Z12 and IPefiarrubia et al.l (120081 ) found that halos 
begin to have their stars stripped after losing 90% of their 
halo mass after infall. 

We use the difference in u max at infall and v max at 
z = for a given VL2 subhalo to estimate the amount of 
mass that it has lost since infall. Following the results of 



IPefiarrubia et al.l (|2010l) . 

W(z = 0) = 2° A x°- 3 (2) 

Vinfall (1 + x) - 3 

for a subhalo with a DM density slope, 7 = —1.0, 
where x = mass(z = 0)/mass(z =infall). A slope of 
-1 is roughly the slope in the centr al regions for all ha- 
los formed in DM-o nly simulations (|Navarro et alJ[T997t 
iSpringel et al.ll2008l ). and we verified that the inner den- 
sity slopes in the VL2 subhalos are also consistent with 
-1. We use the above equation to find (1) all subhalos 
that lose more than 97% of their infall mass (x = 0.03) 
and (2) all subhalos with Ui n f a ii > 30 km s _1 that lose 
more than 90% of their infall mass [x = 0.1) and have 
pericentric passages under 20 kpc. We consider these two 
populations of subhalos to be "destroyed." 

3.2. Identifying Dark Subhalos 

As wc will show below, even after considering destruc- 
tion due to tidal effects on the VL2 subhalos, 80 satellites 
remain at z = from our original sample of 125 with 
« pca k > 20 km s _1 . These 80 subhalos should be bright 
enough to have been detected (see next section), but 80 
is much larger than the number of luminous satellites 
observed in the MW or M31 to date. However, many of 
these satellites are in the halo mass range that is expected 
to be strongly affected by UV heating, and some of these 
subhalos should therefore be inefficient at forming stars 
and will remain dark. In this section, we apply the re- 
sults of lOkamoto et all (|2008l ) to identify the subhalos 
t hat may remain dark. 

lOkamoto et alj (|2008D use hydrodynamical simulations 
to identify the characteristic halo mass, Af c h ar , that re- 
tains 50% of the cosmic baryon fraction, /b ar , as a func- 
tion of z. Their simulations adopt a uniform ionizing 
background that accounts for H and He I reionization 
at z = 9, and Hell reionization at z = 3.5. We have 
converted their M c h ar (z) results into w c h a r(z), adopting a 
WMAP3 cosmology and using an overdensity of 200p cr ;t. 
We assume that if a VL2 subhalo has a u pca k value above 
fch a r(zpe a k) it retains enough baryons to be luminous. 
Note that t>char(z) is a virial quantity, but we do not have 
the virial masses of the VL2 subhalos. Rather, we have 
Vmax values. Most satellites have 1.2 < v max /v v i T < 1.8, 
corresponding to a range of concentration^ 10 < c < 40, 
with a mean ^maxA'char = 1-4 for a h alo with c = 20 
(|Bullock et al.ll200ll : IPrada et alJl20H . 

Figure 2 illustrates the affect of applying this model 
to the VL2 subhalo population. Diamonds show the 
^peak(%>eak) values of the VL2 satellites that survive af- 
te r our tidal stripping co nsiderations. To put the v c hai(z) 
of lOkamoto eTafl p008f) into u max space for comparison 
with the VL2 subhalos, we have multiplied v c har(-z) by 
1.4 to derive the solid line in Figure 2 (i.e., we assume 
"max/«char = 1-4, a typical value for a subhalo). The 
shaded region surrounding the solid line shows the full 
range of 1.2 < v max /v cilal < 1.8. 

Clearly, the number of subhalos above v c har depends 
sensitively on the concentrations (and hence u m ax/^char) 

7 A halo's concentration, c = R v i r /R B , relates the the virial 
radius to a scale radius, R B , where the steep, outer density profile 
of the halo transitions to the less steep, inner profile, e.g., where 
7 ~ -2 for an NFW halo. 
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Fig. 2. — Criteria to select dark subhalos. ^ pea k i s shown for VL2 
subhalos that survive tidal stripping, at the redshift of v poa k. The 
solid line shows 1.4^ c har following Okamoto ct al. (2008), as a func- 
tion of z. The shaded region shows varying 1.2 < !)max/»char < 1-8. 
For our analysis, VL2 subhalos below the solid line are considered 
"dark." 



of individual subhalos. Adopting Umax/^char = 1.4 (1.2) 
leads to 21 (40) surviving, luminous satellites at z = 0. 
Certainly scatter in concentration is expected amongst 
the individual VL2 halos, which we have not accounted 
for. Further, concentrations at these halo masses a re ex- 
pected to decrea s e with increasing z (IBullock et al.l 
Eke et all [20011: iZhao et al.l [20031: iMaccio et alT 
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2008 



Prada et al. 1 120121) . A gain, we do not account for the 



change in concentration with z, but note that the trend 
would allow for lower i> poa k halos to lie above w c har with 
increasing z. The overall determination of whether a par- 
ticular halo is likely to be luminous will depend on the 
halo's formation and evolution. For the purposes of this 
work, we assume that all halos above the solid line in 
Figure 2 (21 subhalos with u m ax/w c har > 1-4) are likely 
to be luminous. 

3.3. Assigning Luminosities 

Finally, we assess whether the surviving satellites are 
luminous by using the t>i n f a ii — M sta r relation from Z123 
to assign stellar masses to the VL2 subhalos, 



= 0.018 



M 

In these simulations, M s tar 



Winfall 



(3) 



km s 

oc Ml ir HGovernato et alj|201^ 



For the stellar mass range probed here (M s t ar >10 Mq), the 
slope of the adopted M a tar-M„j r relation lies between values 
commonly adopted (e.g., [Koposov ct al. 2009; Kravtsoy] 120101 : 



Rashkov ct al. 2012). M sta 
re.E- IKlvpin et al.ll20Tlll . 



because v„ 



scales as M 



1/3 



Fig. 3. — fikpc vs My of the VL2 subhalos. The top panel is the 
direct result from VL2 at z = 0, while the bottom panel shows the 
corrected kinematics based on Z12, along with considerations about 
which subhalos are likely to be observable. Filled red symbols arc 
those satellites that should be observable at z = 0. Subhalos that 
are unlikely to survive due to the tidal effects of a baryonic disk 
are marked by circles with an x through them. Empty symbols are 
subhalos that are likely to be dark. Filled black circles should be 
luminous, and do not experience enough stripping to satisfy our 
destruction criteria, but have lost enough mass that stars should 
be stripped and the luminosities should be considered upper limits. 



We use the tight log(M sta r)— My for the subhalos, 
Iog 10 (-^) = -0.38 - 4.63M y , 



M 



(4) 



in Z12 to further assign V-band magnitudes, My. Z12 
showed that this relation produced simulated satellite 
luminosity functions that were in good agreement with 
the MW and M31. Additionally. iMunshi et al.l (|2012[ ) 
have shown that the simulations that equations 3 and 4 
are drawn from yield an excellent match to the z = 
stellar-to-halo ma ss relation from iMoster et all (|2012| ). 
Most importantly, iGovernato et all (|2012f ) demonstrated 
that these mass-to-luminosity relations produced excel- 
lent agreement with held galaxies in the same luminosity 
range as the satellites we examine here. In other words, 
despit e the fact that past simulations have overproduced 
stars dZolotov et all 120091: iGuo et alj|2010t ^ 



2011 



iBrooks et all l201lt lL^rtaerT T2012: 



Sawala et al 



Moster et al 



2012I ). the simulations used to derive equations 3 and 4 
are in excellent agreement with observed stellar-to-halo 
mass relations. There is no indication that the luminosi- 
ties predicted by these relations are too bright. 

Figure 3 shows the resulting My and wikpc for the VL2 
satellites. The top panel shows the results directly from 
the VL2 catalog at z = 0, while the bottom panel shows 
the results after the destruction, heating, and velocity 
corrections considered in this paper. Filled red data 
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points arc those subhalos that are likely to survive tidal 
effects and retain enough baryons to be luminous, circles 
with 'x' through them are those that are likely to have 
been destroyed by tidal effects in the presence of baryons, 
and empty circles are those that fall below the character- 
istic mass to retain baryons and form stars (see Figure 
2). The filled black circles in Figure 3 identify a popu- 
lation of satellites that have lost more than 90% of their 
mass since infall, but do not meet the destruction criteria 
that we outlined above. iPeharrubia et all (|2008l ) found 
that halos that lose more than 90% of their mass begin to 
strip stars (see also Z12). Hence, the luminosities of the 
black points in Figure 3 should be considered upper lim- 
its, as these halos will likely have had stars stripped and 
be fainter than at infall. We note that My = — 10 corre- 
sponds to Winfaii = 20 km s~ 1 . The few subhalos fainter 
than My = — 10 are those that had u p0 ak > 20 km s _1 
but a Vinfaii < 20 km s _1 . Had we adopted i> p0 ak to assign 
stellar masses, all subhalos in Figure 3 would be brighter 
than My = —10. However, whether we adopt v poa k or 
■finfaii has no effect on the number, masses, or luminosities 
of surviving, luminous subhalos (red circles in Figure 3), 
and does not change our conclusions. We adopt Ui n f a ii to 
assign stellar masses because this is the relation defined 
in Z12. Our focus on subhalos with My < —10 in Figure 
3 is sufficient to make a statement about whether bary- 
onic physics can address the massive missing satellites 
problem. 

Most of the brightest (My < —12) subhalos meet our 
criteria for destruction. We verified that, with the ex- 
ception of one massive satellite accreted at z = 0.5, all 
of these satellites were accreted z > 1.5. The majority 
have Wpeak > 40 km s _1 , and by definition have orbits 
that take them within 20 kpc of the center of the parent 
halo. This is consistent with the idea that early, massive 
satellites contribute to the gro wth of the inner stellar halo 
in Milky Way -mass galaxies (jBullock fc Johnston! 120051 
iZolotov etaLl 120091 ). In DM-only runs that neglect the 
presence of the disk, some of the cuspy inner remnants 
of these early, massive satellites are capable of surviving 
to z = 0. 

Before applying the Z12 correction, the VL2 run con- 
tains more than 20 luminous satellites that have vikpc > 
20 km s -1 , completely inconsistent with the satellite 
popu lation of either the MW or M31 (|Brooks fe Zolotovl 
2012). After applying the Z12 correction and consider- 
ing satellites that are likely to be destroyed by baryonic 
physics or remain dark, only 3 satellites with vik pc > 20 
km s" 1 remain. These 3 satellites are all more luminous 
than Fornax, the MWs brightest dSph, which we discuss 
further in the next section. 

It can be seen from Figure 3 that applying the Z12 cor- 
rection over-corrects a small number of subhalos, result- 
ing in negative velocities. These are a population of halos 
that, by definition, begin with t> max > 20 km s _1 at in- 
fall, but lose so much mass that they have very low «ik pc 
values at z = 0, and result in an overcorrection. While 
Figure 3 already identifies these as halos that have lost 
at least 97% of their mass since infall, wc verified that 
they actually lost more than 99.9% of their mass after 
infall, and that the majority have tidal radii less than 
1 kpc, and thus can safely be associated with destroyed 
subhalos. 



4. DISCUSSION 

Two conclusions can be drawn from Figure 3. The first 
conclusion is that a correction such as that suggested by 
Z12 is required to bring the masses (and hence veloci- 
ties) of the predicted subhalo population in line with ob- 
servational results. Neglecting the baryonic effects that 
reduce the central masses of subhalos will inevitably lead 
to a population of subhalos substantially more massive 
(«ik P c > 20 km s" 1 ) than any of the dSphs observed 
around the Milky Way. The second conclusion is that 
both tidal destruction in the presence of a baryonic disk 
and UV heating must be considered in order to bring the 
total number of luminous satellites in line with observa- 
tions. All subhalos with My brighter than —10 should 
be bright enough to be detected around our MW. The 
uncorrected VL2 catalog would suggest that more than 
100 detectable subhalos should exist, grossly inconsistent 
with the classical dSph population (roughly a dozen in 
the MW and two dozen in M31). Considering destruc- 
tion mechanisms alone only reduces this population to 
80 satellites. An additional correction that assumes sup- 
pression of star formation is necessary to further reduce 
the number of luminous satellites. 

Considering the effects of tidal destruction and heat- 
ing can substantially reduce the number of luminous 
satellites that are predicted to survive at z = around a 
Milky Way-mass galaxy. The exact number depends on 
the assumptions adopted, particularly for the number of 
dark satellites. The model adopted in this work is not 
intended to be conclusive, but rather to motivate more 
rigorous work on this topic. Ideally, future work will 
adopt a semi-analytic model that follows the growth and 
merger history of individual halos to determine if they 
are massive enough to retain baryons and form stars. 
More work is also needed to understand the influences of 
the disk on the survival of substructure. Previous work 
has examined the influence of tidal stripping, though 
not always with the added presence of a galaxy disk 
(iTavlor fc Babulll200l"t IStoehr et alll200l lHavashi et all 



l[Sto 
]J20 

Read et all I2006ri IPeharrubia et all 1200c 



20031: IKravtsov et all 120041: iKazantzidis et alT" 



3b} 

Romano-Dfaz et all [2010bt iNickerson et all 120111) . In 



2004; 



2010; 



this work, we have emphasized the effect of tidal 
stripping on satellites in the presence of a host disk. 
The Z12 correction, however, neglects disk shocking 
that occurs whe n a subhalo passes dir e ctly through the 
baryon ic disk (jTavlor fc Babull 120011: ID'Onghia et~aTl 
l2010al) . and will lead to even faster disruption of a 
subhaloB Tidal he a ting of subhalos should also occur 
dGnedin et all 119991: iMayer et all [20011: ID'On ghia et all 
I2010U IKazantzidis et al.l I2011D. b ut requires very high 
resolution to capture (|Choi et al.l 120091) and is unlikely 



to be accounted for properly in the Z12 correction. 
Hence, all of these processes require more thorough 
study to fully understand the influence of baryons on 
the evolution of satellites. Finally, one galaxy simulation 
alone is not sufficient to understand whether baryonic 
processes can solve the MSP. A statistical sample of 
MW-mass realizations must be used to quantify the 

9 Z12 did have subhalos that experienced disk shocking, but 
these few halos were disrupted so strongly compared to other sub- 
halos with larger pericenterric distances that the correction con- 
siders them outliers. 
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theoretical predictions. 

By the same token, one observed galaxy alone is not 
sufficient to understand whether baryonic processes can 
solve the MSP. Looking at two galaxies, the MW and 
M31, it is already clear that the subhalo population of 
Milky Way-mass galaxies can vary significantly. The sur- 
viving luminous VL2 satellites (shown as red circles in 
Figure 3) include 3 satellites with luminosities and veloc- 
ities larger than the dSph population of the MW. How- 
ever, M31 does have such galaxies, and the surviving 
VL2 subhalos would more closely resemble the luminous 
satellite population of M31. M31 has at least 4 luminous 
satellites with magnitudes brighter than My < —14. 
Three of these arc dwarf ellipticals, a population that 
is missing in the MW. All 3 of the VL2 subhalos were 
accreted at z > 1, have orbital pericenters under 50 
kpc, and are unlikely to retain gas, making them prob- 
ably look very similar to the dwarf ellipticals of M31. 
We also note that many of the VL2 satellites identi- 
fied as luminous at z = reach u pea k at relatively low 
z (see Figure 2), and may have extended star forma- 
tion to quite low z. Z12 also found that satellites with 
My < —8 had extended star format ion histories, consis- 
tent with observational estimates (IGrebel fc Gallagher! 
l200i IDellenbusch etHI [20081: IWeisz et all 120 111) . We 
anticipate that the more recently discovered ultra-faint 
satellites s hould have had the ir star formation truncated 
at high z ([Brown et al.ll2012D . 

In comparing our results for the luminous satellites 
of VL2 to those of the MW and M31, we should con- 
sider the impact of halo host mass on satellite popu- 
lations. It has been shown that the number of sub- 
halos at a given mass scales with host mass, since the 
mass assembly of DM halos is self-similar ([Stewart et al.l 
120081: iFakhouri et~aT] 12010ft . Several works have ar- 
gued that the halo mass of M31 is nearly twice as 
mass i ve as the MW's halo mass (e.g. , iKallivavalil et all 
l200l IGuo et all [20101 IWatkins et all 12010ft . suggesting 
that M31 should have both a brighter satellite popula- 
tion, and more massive satellites, than the MW. While 
the exact halo mass of the MW is still unknown, re- 
cent estimates ( e .g.. [Smith et all [20071: IXue et al.l 120081: 
IGuo et all [2010: IGnedin et al.l 12010ft find a halo mass 
range of 0.7 - 2.0 x 10 12 M Q . The virial mass of VL2 
(1.7 x 10 12 M Q ) is near the upper range of the MW's mass 
estim ates, but perhaps mu ch closer to the virial mass of 
M31 (jWatkins et al-lhoiOft . particularly if M31 is twice 
as massive as the MW. Therefore, the resemblance of the 
VL2 bright satellite population, i.e., its 3 satellites with 
Mv < —14, to the bright satellites of M31 is perhaps in- 
dicative that the halo mass of VL2 more clo sely matches 
the halo mass of M31 than t he MW (see also I Wang et al.l 
I20T3 IVera-Ciro et al.ll20r2ft . 

Because this simulated population is a better match to 
M31 than the MW, we should also consider how com- 
mon the MW satellite distribution is. A few authors 
have recently begun exploring this question. On the 
observational side, searches for satellites in the SDSS 
around galaxies as luminous as the M W find ^10 or fewer 
satellites more luminous than Fornax dBusha et al.ll201 lh 
Guo et all 120111; ILares et al.l I201H IStrigari fc Wechslerl 



los brighter than My < — 14 consistently exist in the six 
MW-mass halos in the high resolution Aquar ius DM-only 
runs. Usin g a larger statistica l samp le, both I Wang et al.l 
((2012ft and lPurcell fc Zentnerl ([20llh find a 10-20% prob- 
ability of finding a subhalo population similar to the 
MW's at a halo mass of 10 12 M Q . This suggests that the 
MW is somewhat rare, though not exceedingly so. The 
rarity of the MW's subhalo luminosity function seems to 
be linked to a gap in luminosities between Fornax (at 
My = -13.4) and the SMC (at M v = -16.8), though 
more work needs to be done to quantify how common 
the gap is, and if having a pair of galaxies as bright as 
the Magellanic Clouds has any influence of the existence 
of such a gap. 

Despite the caveats listed above, it is clear that 
baryons offer a promising solution to solving the MSP, 
without an additional form of warm or self-interacting 
dark matter. Certainly heating and destruction have 
long been considered potential solutions to bring 
the number of predicted satellit es in line wi t h the 
number of obs e rved satellites dBullock et al.l 120001: 
Somervilld [2001 iKravtsov et all [20041: ISimon fc Gehg 
20071: IKoposov et all 120091: INickerson et al.l 120111 
Penarrubia et alj 120121 ). The key point, however, is 
that even if the number of luminous satellites can be 
brought into line, the overall distribution in masses 
(and obse rved velocities) of those satellite s is usually 
too high (jBovlan-Kolchin et al.l [20TT[ [20121) . The Z12 
model builds on this previous work by considering 
the additional impact of baryons on the dark matter 
density slopes of satellites, finding shallow profiles 
in satellites brighter than My < —12, which further 
reduces their mass. When combined with more effective 
tides after infall dTavlor fc Babull [200lt IStoehr et al] 
[2001 IRead et al.ll2006bt IPenarrubia et alj 12010ft . both 
the number and masses of the subhalo population can 
be brought into agreement with observations for the 
first time. 

In fact, the model utilized by Z12 goes beyond solv- 
ing the MSP or TBTF problems. The model in Z12 has 
also been shown to creat e bulgeless dwarf disk galax- 
ies ([Governato et al.ll2010ft and s maller bulges in higher 
mass galaxies jBrook et al.ll2012l ). This is because winds 
driven by supernovae preferentially remove low angu- 
lar mom entum material fr om galaxies (Christensen et 
al. 2012: IBrook et alj I201T1 ). Importantly, these winds 
are naturally driven (rather than artificially inserted) 
when supernovae deposit thermal energy into a high den- 
sity medium, and create highly overpressurized bubbles. 
Hence, a model must be able to reproduce the high den- 
sities found in star forming molecular clouds (> 10 — 100 
atoms/cc) to naturally drive such winds. Critically, the 
expanding overpressurized bubbles are able to drive rapid 
fluctuations in the potential well of a galaxy. The cumu- 
lative effect is to expand the central orbits of dark matter, 
transforming an initially steep, cuspy DM density pro- 
file int o a shallow, cored profile (e.g.. IRead fc Gilmord 
2005 1 : IMashchenko et all [20061: IPontzen fc Goyernatc 
20121 IMaccio et all 120121: IPontzen fc Governatol 120121 : 



2012ft . Adopting a semi-analyt ic model to assign stellar 



Tevssier et all 12012ft IGovernato et al.l (|2012l) and 



Oh et all ([201 lft showed that the density profiles of sim- 



masses, IVera-Ciro et alj ([2012ft find that 2 to 5 subha 



ulatcd galaxies within this model are in excellent agree- 
ment with the high resolution observations of field galax- 
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In summary, the model adopted by Z12 not only yields 
a population of satellite galaxies in agreement with MW 
or M31 observations, potentially solving both the "miss- 
ing satellites" and "too big to fail" problems, but it 
also simultaneously addresses other small scale problems 
within CDM. Hence, a unified baryonic solution remains 
viable to solve the small scale crisis of CDM. Future work 
must quantify the impact of baryons on the DM struc- 
ture of galaxies, and make testable predictions across all 
galaxy masses. 

5. CONCLUSIONS 

It has long been appreciated that the observed pop- 
ulation of Milky Way satellites is at odds with the dis- 
tribution predicted by dark matter-only simulations. In 
this article, we have argued that effects associated with 
baryonic physics can reconcile the results of such sim- 
ulations with observations, without the need for dark 
matter that is warm, self-interacting, or with properties 
that are otherwise different from those of the standard 
cold and collisionsless paradigm. In particular, super- 
nova feedback in luminous satellites can reduce the den- 
sity of dark matter in the inner volumes of these systems, 
while the presence of a baryonic disk can enhance the de- 
gree of tidal stripping that takes place. The combination 
of these effects leads to a reduction of the masses of the 
predicted subhalo population, bringing the overall num- 
ber of satellites into concordance with that observed. In 
particular, in the Via Lactea II sample considered, we 
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found that these effects reduced the number of subhalos 
with masses greater than 10 s M Q from more than 20 to 
only a few. 

After determining the distribution of massive dark 
matter subhalos predicted to be present in a Milky Way- 
mass galaxy, we turned our attention on the question 
of how many of these objects are likely to host lumi- 
nous satellites. For reasonable assumptions regarding 
the stellar-to-halo mass relationship and for the criteria 
for destruction via tidal stripping, we predict a luminous 
satellites population that is in adequate agreement with 
both the Milky Way and M31. 

While the work presented here is not intended to rep- 
resent the final word on this topic, we have shown that 
baryonic effects can lead to a population of satellites 
around Milky Way-mass galaxies that is is good agree- 
ment with observations. This strongly reduces the moti- 
vation for warm or self-interacting dark matter scenarios 
and represents yet another success for cold dark matter. 
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